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ABSTRACT 

We report the discovery and characterisation of a deeply eclipsing AM CVn-system, 
Gaial4aae (= ASSASN-14cn). Gaial4aae was identified inde pend ently by the All-Sky Au¬ 
tomated Survey for Supernovae (ASAS-SN; IShappee et al. i 120141) and by the Gaia Sci¬ 
ence Alerts proje ct, during two separa t e outbursts. A third outburst is seen in a rchival Pan- 
STARRS-1 (PS 1: ISchlaflv et al. Il2012t fTonrv et al. ll20H iMagnier et al. 11201 3l) and ASAS- 
SN data. Spectroscopy reveals a hot, hydrogen-deficient spectrum with clear double-peaked 
emission lines, consistent with an accreting double degenerate classification. We use follow¬ 
up photometry to constrain the orbital parameters of the system. We find an orbital period 
of 49.71 min, which places Gaial4aae at the long period extremum of the outbursting AM 
CVn period distribution. Gaial4aae is dominated by the light from its accreting white dwarf. 
Assuming an orbital inclination of 90° for the binary system, the contact phases of the white 
dwarf lead to lower limits of 0.78 M 0 and 0.015 M 0 on the masses of the accretor and donor 
respectively and a lower limit on the mass ratio of 0.019. Gaial4aae is only the third eclipsing 
AM CVn star known, and the first in which the WD is totally eclipsed. Using a helium WD 
model, we estimate the accretor’s effective temperature to be 12900 ± 200 K. The three out¬ 
burst events occurred within 4 months of each other, while no other outburs t activit y is seen 
in the previous 8 years of Catalina Real-time Transient Survey (CRTS: [Drake et aT71l2009b . 
Pan-STARRS-1 and ASAS-SN data. This suggests that these events might be rebrightenings 
of the first outburst rather than individual events. 

Key words: Binaries: eclipsing. Stars: Cataclysmic Variables 


1 INTRODUCTION 


AM Canum Venaticorum (AM CVn) stars are a rare class of com¬ 
pact hydrogen-deficient interacting binaries, comprised of white 
dwarfs (WDs) accreting He-rich materi al from low mass degener¬ 
ate or semi-degenerate companions (see lNeleman s 2005; S olheiml 
for recent reviews). The orbital periods of these systems 
range from 5 to 65 min. This implies highly evolved compo¬ 
nents and makes them, along with their ultra-compact X-ray bi¬ 
nary equivalents, one of the most compact classes of binary sys¬ 
tem known. The prototype system for the class of object was 
discovered in 1967, an d has an orbital period of 17 min ISmakl 
1 19671 : IPaczvnski Ifl967l) . Since then, 43 confirmed AM CVn sys¬ 
tems have been discovered dLevitan et al, l2015h . Not only are these 
systems interesting as on e of th e pos sible end points for binary 
WD evolution iNelemans et al. 1120011) . they are also potentially 
strong sources of gravitational wave emission due to their com¬ 
pact configurations dNelemans 120031). and they may be the progen¬ 
itors of peculiar “dot la” supernovae i Solheim & Yungelson 120051: 


iBildsten et al. l2007l:llnserra et al. Il2015t). 

As binaries, AM CVn systems can yield detailed informa¬ 
tion on the masses and radii of the two components if eclipses 
and radial velocity variations can be observed. Eclipsing systems 
in particular offer the possibility of measuring full system parame¬ 
ters, including inclination and component masses, from time-series 
photometry alone. The most robust results come from systems in 
which the white dwarf is totally eclipsed. The extreme mass ra¬ 
tios of AM CVns mean that the likelihood of observing such sys¬ 
tems is low and currently only two eclipsing AM CVn systems 
are known. SDSSJ0926+3624 was the fi rst eclipsing AM CV n star 
to be discoveredj Anderson et al. Il2005l : ICopperwheat et al~ll201 ll : 


ISzvprvt et al. 112014 ), however its WD is only partia lly eclipsed. A 

second partial eclipser (PTF1 J191905.19+481506.2, iLevitan et al. I 
Emi) was recently discovered, but it only eclipses the edge of the 
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disc and not the WD, and so cannot be used for parameter determi¬ 
nation. 

Determining the nature of the secondary (donor) star is criti¬ 
cal to our understanding of the past evolution of the system, since 
the three binary evolution channels proposed to form AM CVn 
stars are best distinguished by the state of the donor star at the 
onset of mass transfer. If the primary WD is accreting from an¬ 
other He-rich WD, the binary must have undergone two common 
envelope events in the past to reduce it to the observed compact 
configuration. This is known as the d ouble degenerate channel 
dPaczvnski 1 1967l : lFaulkner et al. 1 1 9721) . Alternatively, if the donor 
is not fully degenerate at the time when it leaves the second com¬ 
mon envelope, the donor will be more massive than in the case 
of the double-degenerate channel. Mass loss will cause i t to be¬ 
come increasingly degenerate as the binary evolves dSavoniie et al. I 
Il986l : llben & Tutukov II 19871) . At the longest observed orbital peri¬ 
ods (i.e., the oldest AM CVn systems), the two channels become in¬ 
distinguishable. The donor is predicted to reach the same near-zero 
temperature, low entropy configuration in both cases dDelove et al. 1 
120071) . A third possibility is that the binary may start mass trans¬ 
fer as a hydrogen-rich cataclysmic variable. Such a system could 
evolve to an AM CVn star if the donor star had already started to 
evolve by the time mass transfer starts, and results in a hotter, more 
massive donor star and traces of hy drogen may be expected in such 
systems (Podsiadlowski et al. 2 0031) . 

In this paper, we present follow-up observations and pre¬ 
liminary modelling of the AM CVn system Gaia l4aae (RA = 
16:11:33.97, Dec= +63:08:31.8. iRixon et al. 1 120141) . Gaial4aae 
was first detected in outburst by the All -S ky Au tomated Survey 
for Supemovae (ASAS-SN. Ishappee et al. II2OI4I ) at V = 13.6 on 
2014 June 14, who gave it the designation ASASSN-14cn. This 
was before the formal start of the Gaia Science Alerts project. 
However, about 2 months later, Gaial4aae underwent a second out¬ 
burst, which was detected by Gaia on 2014 August 11 at G = 16.04, 
during the science commissioning phase. As this was significantly 
(1.52 mag) brighter than the historic Gaia magnitude of the source 
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at this position, it was identified and announced as a Gaia science 
alerQ. 

Gaial4aae was discovered as part of the Gaia Science Alert s 
(GSA) project (Hodgkin et al. in prep.; Wyrzyk owski et al~ll2012h . 
which aims to identify such photometric transients in the Gaia 
satellite data, and publicly announce their discovery on a rapid 
timescale. Gaia is scanning the entire sky at sub-milliarcsecond 
resolution with precise photometry and astrometry down to a lim¬ 
iting magnitude of G~20 (G is the Gaia white light bandpass; 
[jordi et aklboioh . Over the five year mission each position on the 
sky will be observed on average 70 times. These repeated obser¬ 
vations of the entire sky mean that alongside the primary science 
mission of Gaia, to provide spatial, kinematic and physical param¬ 
eters for a billion stars in the Milky Way, the satellite will also ob¬ 
serve many transient and time-domain phenomena, which will be 
explored systematically by the GSA project. 

Many of the known AM CVn systems d isplay outbursts and 
super-outbursts in th eir long-term light curves J Levi tan et al. 20 T 5 I: 
iRamsav et al. Il20i~2h . during which they brighten by 3-4 mag over 
timescales of 1-2 days and last weeks to months. Currently, it ap¬ 
pears that about 60 per cent (27/44) of the known AM CVn systems 
display outbursts iLevitan et al. [|2015ll . Gaia will play an impor¬ 
tant role in the discovery of new cataclysmic variables (CVs), both 
in outburst, and also through their decrease in magnitude during 
eclipses. From pre-launch simulations, we expect that ~1000 new 
CVs, including a number of evolved systems and AM CVn systems 
will be found by Gaia over its mission lifetime. 


2 OBSERVATIONAL DATA 

A 300s long-slit spectrum of Gaial4aae was taken on the night of 
2014 October 13 (MJD 56943.88751 at the mid-point of the expo¬ 
sure), when the system had returned to its quiescent state (i = 18.74 
± 0.02 mag). This spectrum was obtained using the Auxiliary-port 
camera (ACAM), with the V400 grating, on the 4.2-m William Her- 
schel Telescope (WHT). The data were reduced within IRAF in the 
standard fashion. The extracted and calibrated spectrum has a reso¬ 
lution of ~ 12 A and a S/N of ~20 in the continuum. The spectrum, 
plotted in Fig. Q] shows clear He emission lines, but no detectable 
H lines. The emission lines are broad (Full width at half maximum 
(FWHM) = 2415 ± 100 km/s) and display double-peaked profiles, 
which are typical of AM CVn stars, revealing the presence of an 
He-dominated accretion disc. Based on this spectrum, we classi¬ 
fied Gaial4aae as an AM CVn system. The spectral classification 
and similarities to other eclipses motivated the acquisition of fur¬ 
ther photometry. The peak velocities of the two emission compo¬ 
nents are at 800 ± 50 km/s relative to rest frame, averaged over 
all detected emission lines. Measurements of the peak separation 
of the individual double peaked lines all agree within 3<r of the av¬ 
erage value. We do not see a sharp central spike between the lines, 
which is observed in many A M CVn stars and thought to origi¬ 
nate on the surface of the WD llMarsh II 19991 ; iMorales-Rueda et al. I 
120031 ; iRoelofs et al. 1 120071 , l2009t) . This might be due to the low 
resolution of the spectrum, although it could also be because of 
the high inclination of this system, as appears to be the case with 
SDSSJ0926+3624 JConnerwheat et afluOl ll) . On the ephemeris 
given in Section[3] the WHT spectrum was taken away from eclipse 
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Table 1. Log of photometric observations of Gaial4aae used in this work. 


Observatory 

Obs. date (UT) 

Filter 

Exposures (s) 

Gaia 

2014 08 11 

G 

45 

ASAS-SN 

2012-2015 

V 

129x180 

Loiano 1.5m Cassini 

2014 10 24 

g 

3x300,91x30 

Telescope + BFOSC 

2014 10 25 

g 

135x30 

Bialkow 0.6m, Poland 

2014 10 18 

BV 

30x120 


2014 10 19 

BV 

37x120 

CIECEM 0.35m, Spain 

2014 10 21 to 
2014 11 18 

clear 

40x180, 8x150 
111x120, 399x90 

pt5m. La Palma 

2014 10 25 

V 

61x60 


2014 10 22 

V 

36x60, 21x120 

0.6m ASV, Serbia 

2014 10 21 

BVRI 

6x300 

Belogradchik AO 0.6m, 
Bulgaria 

2014 10 21 

BVR 

2x 300 

Asiago 1.82m Copemico 

2014 12 11 

r 

169x20 


2014 12 12 

g 

169x20 

4.2m WHT+ACAM 

2014 12 18 

V 

491x5 

Mercator 

2015 01 15 

g r+i 

232x30 

Catalina (historic) 

2005 - 2014 

clear 

107x30 

Pan-STARRSl (historic) 

2010 - 2014 

grizy 

66x30 


at phase 0.22. Thus, the absence of the emission spike is not due to 
the white dwarf eclipse. 

The historic optical and infrared fluxes of Gaial4aae in pre¬ 
sumed quiescence are also shown in Fig. [T] The optical fluxes are 
from the Sloan Digital Sky Survey (SDSS) DR 10 iAihara et al. I 
Imj]), while the infrar ed fluxes are from forced photometry at the 
SDSS source loc ation (Lang 2014) on Wide-Field Infrared Survey 
Explorer (VP/SEt lWright et al. Il2010h images. Ultraviolet (UV) im¬ 
ages ar e also a vailable from Galaxy Evolution Explorer ( GALEX; 
Martin et al . 12005 ) DR 6. The GALEX archive contains three pairs 
of FUV & NUV observations for Gaial4aae, one obtained on 2005 
March 9, and the other two on 2007 May 24. All three observa¬ 
tions had short exposure times, 143-195 s. One of the May 2007 
observations shows the system at a significantly fainter level than 
the other two. While our current ephemeris is not sufficiently ac¬ 
curate to establish the orbital phases of the GALEX observations, 
it is most likely that the system was caught close to the eclipse of 
the primary. The eclipse duration, discussed in Section 3, is Ills, 
shorter than but comparable to the GALEX observations. All fluxes 
have been corrected for extinction towards the source, E(B — V) 
= 0.018 dSchlegel et al. Ill998h . The absolute flux calibration of the 
WHT+ACAM spectrum has been scaled to match the SDSS r and 
i band magnitudes for Gaial4aae. 

The initial determination that Gaial4aae was eclipsing 
was made by the “Centre for Backyard Astrophysics” project 
I lSkillman & Patterson 1 119931 ); who established a preliminary pe¬ 
riod for Gaial4aae of 49.7 min t ide MigueTll20l4 ). Following this, 
an intensive photometric monitoring campaign was undertaken for 
Gaial4aae at a number of telescopes, as detailed in Table Q] In 
addition to this, we searched the databases of the Catalina Real¬ 
time Transient Survey (CRTS; Drake et al. ||2009h, Pan-STARRS- 1 


time Transient Survey (CR T S; Drake et al. 2UU9 ), Fan-STAR RIS-1 

(PSLlMagnieretal. 12013; Schlaflv et al. l2012l:lTonrv et al. l2012h 


and ASAS-SN (Shappee et al. 20141) for pre-discovery images 
covering the position of Gaial4aae. The cadence of the CRTS data 
is relatively low, but during those observations no outbursts were 
observed. The average quiescent magnitude in CRTS for Gaial4aae 
is 18.64 ± 0.14 mag. PS1 detected an outburst of Gaial4aae on 
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Wavelength [A] 

Figure 1. WHT+ACAM spectrum of Gaial4aae taken on 2014 October 13 during quiescence, showing double-peaked He emission and an absence of H lines. 
The historical GALEX and SDSS photometry are also shown as red points; the fainter GALEX magnitudes probably cover an eclipse. The blue and magenta 
lines are T’eff = 12700 K and 13100 K respectively He-atmosphere models fitted to the optical flux and two different epochs of UV flux. The top-right inset 
shows the spectral energy distribution (SED) fit with the WISE data included, the bottom-light inset shows a zoom in of the He I A6678 line in velocity space. 


2014 July 7, when it reached 15.38 mag in /-band, compared to 
18.74 ± 0.02 mag in quiescence. Two eclipses of Gaial4aae are 
also visible in the PS1 data. ASAS-SN has many upper limits for 
the light curve and detected the decline of the outburst they dis¬ 
covered, as well as some data on the second outburst, but only 
place limits on the third outburst. The combined light curve for 
Gaial4aae spanning 8 years of PS1, CRTS, ASAS-SN and Gaia 
data is shown in Fig. [2] The first Gaia data point shown is the av¬ 
erage of the 1.5 days of data Gaia had observed before the outburst 
was discovered. This may already included some of the rise of the 
outburst and thus be higher than the true historic magnitude. 

From the combined light curve, it appears that Gaial4aae un¬ 
derwent at least three outbursts between 2014 June and September. 
The first outburst was seen by ASAS-SN on 2014 June 14. The sec¬ 
ond outburst was seen by PS1 in /-band on 2014 July 7 and ASAS- 
SN on 2014 July 8. The limits measured by ASAS-SN between 
2014 June 20 (6 days after the first outburst) and June 27 (9 days 
before the second outburst) rule out the possibility that the first and 
second outburst are in fact one continuing event. The third outburst 
of the system was caught by Gaia on 2014 August 13, and is con¬ 
strained by the Gaia historic data 1.5 days prior, the PS1 detections 
of the system in quiescence in i band 24 days prior and in z band 7 
days after, as well as ASAS-SN limits 1 day after, suggesting this 
outburst had a short duration. 

A number of follow-up studies were conducted. Imaging ob¬ 
tained with the pt5m. La Palma (Hardy et al. in prep.) was re¬ 
duced using the ULTRACAM pipeline {iDhillon et al. l2007h . while 
for all other instruments with the exception of WHT+ACAM, the 
data were debiased and flatfielded using standard techniques. The 
ACAM data was taken using a small CCD window and a fast 
readout mode; as no suitable Hatfield or bias frames were avail¬ 


able, these calibrations have not been applied. However, as we are 
performing differential photometry over a small area on a single 
night, this should not affect our results significantly. ASTROM- 
ETRY.NET dLang et al. Ihoid ) was run on each image, excluding 
the Cassini+BFOSC, Asiago 1.82m Copernico and WHT+ACAM 
data, to re gister it to a common Wo rld Coordinate System. SEX- 
TRACTOR jBertin & Arnouts Ill996h was used to detect, deblend 
and measure the instrumental magnitudes of all sources in the 
field. Finally, the list of sources detected in each image was up- 
loaded to the Camb ridge Photometry Calibration Server (CPCSq 
iWvrzvkowski et al. 1120130 . which calibrates all the data from dif¬ 
ferent telescopes to a common photometric system. To measure the 
magnitude of Gaial4aae on the BFOSC, Asiago and ACAM im¬ 
ages, we used c o-loc atedJist-driven differential photometry as de¬ 
scribed in llrwin et al. 1 ( t2007h . using the CASUTools package, 
yielding a precision of 15-18 millimag for Gaial4aae while out 
of eclipse. The comparison stars were checked and found to be 
photometrically stable. To correct for light travel times, we con¬ 
verted the MJD (UTC) times of all data to the barycentric dynami¬ 
cal timescale (TDB). 


3 ANALYSIS 

In order to estimate the WD temperature, we assume that the con¬ 
tribution of the accretion disc to the GALEX FUV and NUV fluxes 
is negligible, and fit the three ultrav iolet observations with helium- 
atmosphere models from lKoesterl ( bold) , as shown in Fig. [T] We 


2 gsaweb.ast.cam.ac.uk/followup/ 
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Figure 2. Historic light curve from CRTS (black crosses) and Pan-STARRSl (filled circles) spanning eight years of observations. Pan-STARRSl clearly 
detected outbursts from Gaial4aae in 2014, and appears to have seen two eclipses. The ASAS-SN and Gaia detections are shown by the turquoise diamonds 
(and limits as grey triangles) and a purple stars respectively. This historic light curve begins on 2006 January 5 and ends on 2015 Marchl5. 


estimate the contribution of the accretion components from the r- 
band light curve (discusses below), which is consistent with the as¬ 
sumptions used in the DB model. The two sets of “bright” GALEX 
FUV and NUV fluxes are consistent with effective temperature es¬ 
timates of T e ff = 12900 ±200 K (from T e g = 12700 K, magenta 
line and 13100 K, blue line, respectively). DB white dwarfs have 
very weak lines at such low temperatures, and thus are not de¬ 
tectable given the much stronger emission lines at this resolution, 
which might explain the lack of broad WD absorption features in 
the spectrum. Adopting a primary mass of Mi = 0.78 Mg, cor¬ 
responding to the lower limit from the light curve fit (see below), 
implies a radius of Ri = 7.4 4 x 10 s cm (using the cooling models of 
iHolberg & Bergeron II2006I) . and hence a distance of 225 ± 10 pc. 
There appears to be an IR excess in the WISE photometry when 
compared to the He-atmosphere model. The IR excess is unlikely 
to be due to outbursts as the WISE photometry is from observations 
taken over 2 weeks separated by 6 months. The first set of WISE 
data was taken over a period from 2010 July 17-23, while the sec¬ 
ond set was taken over 2010 December 23-29; there are 7 CRTS 
measurements during the first set of WISE observations which con¬ 
strain the system to be in quiescence. The cause of the WISE flux 
excess is unclear. 

The ephemer is of Gaial4aae was first d etermined by fitting a 
light curve model JCopperwheat et ah 2m3) to all the photometric 
data divided into 16 night-long chunks. The model is composed of 


a WD, accretion disc and a bright-spot where the gas stream hits 
the disc. The model took into account the finite exposure lengths 
of the images, including their readout time, by computing over 
sub-steps in each exposure. We found the ephemeris of Gaial4aae 
to be 

BMJD (TDB) = 56980.0557197 (13) + 0.034519487 (16) E, 

where the zero phase corresponds to the mid-point of the eclipse, 
based on the time series data from Loiano, Asiago, and WHT. The 
time of zero phase was chosen to give minimal correlation between 
the two fitted parameters and the quoted uncertainties are the lcr 
errors. At present the estimate of the ephemeris suffers from a few 
caveats. First, the long (30s or 20s) exposures used for the Loiano 
and Asiago data, and the small number of eclipse times used (3, 1 
and 2 from Loiano, Asiago and WHT, respectively) are not ideal. 
Second, none of these instruments are built for precise timing and 
may suffer from systematic s. 

We also modelled the high cadence ACAM light curve shown 
in Fig. [3] The light from this system is dominated by the WD, with 
a small contribution from the disc and bright spot. It is estimated 
that the accretion components (bright spot plus accretion disk) con¬ 
tribute ~30 per cent of the r-band flux, although of course this 
component is variable. The pre-eclipse “hump” which originates 
in the bright spot seems unusually variable, and sometimes can 
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Figure 3. Top: Observed r-band WHT+ACAM light curve for Gaial4aae (points) with the best fitting model (lines) comprising of a WD (which is the main 
contributor to the light), accretion disc and bright-spot where the gas stream hits the disc. Bottom: Zoom in around the echpses of the light curve shown above. 


barely be seen, although this may be due to severe flickering. These 
are aperiodic brightn ess var iations with_characteristic timescales of 
seconds to minutes jMiddleditch & Cordova Il982l) . The amplitude 
of the flickering exceeds the noise and limits the current model fit. 
Further observations are required to average the flickering out. The 
eclipses are sharp-sided and deep, and the mid eclipse depths reach 
around 2 mag. In order to constrain the scaled white dwarf radius, 
7T = Ri/a (where Ri is radius of the primary and a is the bi¬ 
nary separation), we determined the phase of the WD eclipse to be 


AT> = 0.0373 ± 0.0005 from our model fit. The ingress and egress 
phases were deduced from the parameterised model of the binary 
fitted to the WHT+ACAM light curve. This gives us n as a func¬ 
tion of the mass ratio q and the inclination i. If we then assume a 
WD mass - radius relation, we can solve for Mi and M 2 using q, 
ri and the orbital period using Kepl er’s laws. Here, we assume th e 
relation of P. Eggleton as quoted in lVerbunt & Rappaport I d 1988h . 
scaling the relation by a factor of 1.05 to account for the finite tem¬ 
perature of the WD. 
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There are a range of parameters which fit the current data with 
our model. The model fits shown in Fig. [3] are for i = 88°. The 
lower limits to both Mi and M2 correspond to i = 90°, q = 0.019, 
n = 0.026 and a = O.413R 0 . The lower limit on Mi ~ 0.782 Mg is 
consistent withtiiejiverage mass of WDs in cataclysmic variables 
dZorotovic et al. Il201ll) . For i = 90°, the companion star M 2 has a 
mass of 0.015 M@, which is consistent with e xpectations for a near¬ 
zero entropy donor at a period of ~50 min ( jPelove et al. 112007 ). 
For a lower inclination model with i = 80°, q = 0.133, n = 0.013, 
a = 0.488 R 0 , Mi and M2, increase to 1.159 M@ and 0.154 Mq 
respectively. From our current data, we are unable to derive a se¬ 
cure value of the mass ratio q, due to the flickering and the weak 
bright-spot. Flence, we cannot select between the low mass, highly 
degenerate donor stars characteristic of the double WD route as 
found for i ~ 90°, and more massive hot donors that one might ex¬ 
pect from the post-CV route (i ~ 80°). Future high cadence, high 
S/N observations over multiple orbits might allow us to measure 
the bright-spot features and break the degeneracy in our derived 
parameters. 


4 DISCUSSION 


The orbital periods of AM CVn stars are thought to increase as 
mass is transferred from donor to accretor, leadin g to a decreasejn 
the ra t e of mass transfer as th e system evolves dTsugawa & Osakil 
119971 : Ne lemans et aTTIboOlh . Thermal instabilities are expected 
and often observed in AM CVn He accretion discs with interme¬ 
diate mass-transfer rates, and these are sometimes seen as dwarf 
nova (DN) type outbursts (Tsugawa & Osaki 1l 1 9971) . Intermediate 
mass-transfer rates are though t to occur for systems with orbital 
periods of 20 min to ^ 40 min dRamsav et al. |[2()T~3 : lLevitan et al. I 
120151 : iNelemans ll2005h . Longer period objects (P or h > 40 min) 
are thought to have low mass transfer rates and stable cool discs, so 
that these should not have ou tburst s, which is mostly confirmed by 
observations dRamsav et al. II 20 I 2 I) . However, the low mass trans¬ 
fer rate could also mean that the intervals between outbursts are 
very long, so we ha v e simply not obser v ed that many outbursts 
(Levitan et al. 1 1201 5l : iKotko et al, 1 120121 : ICannizzo & Nelemans I 
l2015t) . 


Interestingly, Gaial4aae has experienced three outbursts 
within only three months, while no outbursts were detected in 
~8 years, although we cannot rule out that some could have oc¬ 
curred during gaps in data coverage. Thus to see three outbursts 
in just a few months, suggests that they are likely to be “rebright¬ 
ening” outbursts (also known as echo outbursts), rather than in¬ 
dependent events. Multiple rebrightenings are frequently observed 


(e.g. Patterson et al. 1998|; 

Shears et al. 20121; Kato et al. 20141; 

Meyer & Mever-Hofmeister 

2015). Echo outbursts are very simi- 


lar to “normal” dwarf nova outbursts, except that they happen in 
quick succession in a system with otherwise few observed out¬ 
bursts, and they always happen on the decline from a superout¬ 
burst. From Fig.[2]it can be seen that each outburst reaches a lower 
peak magnitude than the previous outburst, consistent with echo 
outbursts, where overall, the target is fading, but it has a few echo 
outbursts following the superoutburst. In between the rebrighten¬ 
ings it fades to near-quiescence. WZ Sge stars and the outbursting 
AM CVn stars, such as Gaial4aae, both have low mass transfer 
rates and extreme mass ratios, which ar e likely to i mpact on the 
duration and frequency of outbursts. iLevitan et al. I d2015h inves¬ 
tigates the correlation between orbital period and outburst recur¬ 


rence time, by extrapolating to rare, long outbursts for long period 
systems. For our system, with a period of 49.7 mins, they predict 
outbursts to recur every ~10 years, although this does not consider 
rebrightenings. 

It is somewhat surprising that Gaial4aae shows outbursts 
at all, because a system with suc h a lon g orbital period is ex¬ 
pected to have a stable, cool disc l|Solhei mll201(]h . However, re¬ 
cent s tudies by ICannizzo & Nelemans d2015h and I Kotko et al. I 
J20 1 2 I 1 used a disc instability model a nd the observ e d out - 
burst properties of systems, as compiled by iLevitan et al. 1 J2015l l. 
to find that systems with higher mass accretors, have lower 
outburst thresholds and are more likely to undergo outbursts. 
Along with other long period AM CVn stars whic h exper ience 
outbursts, SDSS J090221.35+381941.9 dKato et al. 1 12014|) and 


1.9 (E 

Som 


CSSJ045019.7+093113 JWoudt et al. 12013 D, these authors suggest 
that the transition to a stable disc may happen at longer orbital pe¬ 
riods in some cases (or perhaps not at all). 

The temperature of the WD implies an accretion rate, if 
accretion-heated, of 7 — 8 x 10 -11 M 0 yr -1 for O.75M 0 
iTownslev & Gansicke 1120(7}) . Combined with the masses we de- 
rive , this accre tion rate is more consistent with a degenerate donor 
dDelove et al. i l200~7h . suggesting that the system may have de¬ 
scended from a merging double WD, and that it may have had a 
much shorter orbital period in the past (<10 min). 

For comparison, we can compare the accretion rate implied 
from the WD temperature t o the s tability criteria of the disc insta¬ 
bility model in lKotko et al. lf2oT3) . A disc will be stable in the high 
state if it is too hot and it will be stable in the low state if it is too 
cold. For the system to be unstable, the accretion rate in the disk 
must be between the limits for the critical accretion rate for hot 
(M+ r ) and cold (M~ r ) stable equilibrium accretion rates. For an in¬ 
clination of 90°, with no hydrogen, 98 per cent helium and 2 per 
cent metals, the upper critical rate is M+ r = 5.2xl0~ 9 M 0 yr” 1 

4.3x 10 -12 M 0 yr -1 . The WD 


M 0 yr 1 is 


and lower critical rate is M c 

temperature inferred accretion rate of 7 — 8 x 10 -11 
between the these limits. In fact, for any plausible parameters of 
the disk instability model, the inferred accretion rate is orders of 
magnitude below the hot, stable state, and a factor ~20 above the 
cool, stable state. Thus, Gaial4aae is consistent with the disc insta¬ 
bility models for AM CVn stars, since the accretion rate inferred 
from the WD temperature lies in the unstable regime at this orbital 
period and mass. Our estimate of the accretion rate could be too 
high because the WD temperature at the long period of Gaial4aae 
may be set by simple WD cooling dBildsten et al. 1 20061 )- However, 
it would need to be a factor of 20 lower than we estimate to have 
an accretion rate below the lower critical rate M~. 


In the future, there is a variety of data which will be essen¬ 
tial for fully characterising Gaial4aae. Firstly, more precise, high 
cadence photometry can be used to average out flickering, which 
is limiting the analysis of the light curve at present. It is vital to 
observe the bright spot in the system, as this will precisely pin 
down the orientation, thus allowing the system parameters to be 
accurately calculated. Further spectra may allow us to detect the 
same narrow spikes between the double peaked emission lines that 
are seen in other AM CVn stars. Combined with the phase from 
eclipses, this could allow a definitive proof that the spike originates 
on the accreting WD. Spectra will also provide information on the 
elements present in the system, useful for understanding the evolu¬ 
tionary history of Gaial4aae. Gaia will also provide parallax and 
proper motion, and having an accurate distance to the system will 
allow system parameters, such as the WD temperature, to be bet¬ 
ter constrained. Finally, long term precision timing will be needed 
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to detect the expected period change due to gravitational radiation- 
driven mass transfer. This effect should cause a progressive delay 
in the arrival time of eclipses, but it may be at least a decade before 
this can be detected. 


5 CONCLUSIONS 

Gaial4aae was found as a transient in Gaia data on 2014 August 
11. We undertook spectroscopic and photometric follow-up and 
identify it as an AM CVn system. Gaial4aae is a deeply eclips¬ 
ing system, with the accreting white dwarf being totally eclipsed 
on a period of 0.034519 days (49.71 min). It is the third eclipsing 
AM CVn known, the second in which the white dwarf is eclipsed, 
and the first in which the white dwarf is totally eclipsed. We de¬ 
tected three outbursts over ~4 months. The orbital period places 
Gaial4aae at the long period extremum of the outbursting region 
of the AM CVn distribution. A helium WD model was used to es¬ 
timate an effective temperature of ~ 12900 ± 200 K for the white 
dwarf. We used the contact phases of the WD eclipse to place lower 
limits of 0.78 Mq and 0.015 Mq on the masses of the accretor and 
donor respectively, which correspond to an inclination of 90°, a 
mass ratio of 0.019 and an orbital separation of 0.41 Rq. The deep 
eclipses shown by Gaial4aae, suggest that future observation have 
the potential to lead to the most precise parameter determinations 
of any AM CVn star discovered to date. 
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